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Abstract. Using data from the ER0S2 microlensing survey, we report the discovery of two blue objects with 
irregular photometric behaviour of AV ~ 0.1-0.4 mag on time scales of 20 to 200 days. They show a bluer when 
fainter behaviour. Subsequent spectra taken with the ESO 3.6m telescope show spectral type B4eIII and B2eIV-V 
with strong Ha emission. These objects resemble the Herbig AeBe but also classical Be stars. At this stage, it is 
not possible to distinguish unambiguously between pre-main sequence and classical Be nature. If we favour the 
pre-main sequence interpretation, they are more luminous than the luminosity upper limit for Galactic HAeBe 
stars. The same was found for the HAeBe candidates in the LMC. This might be due to a shorter accretion time 
scale (r — M^/ M), or the smaller dust content during the pre-main sequence evolution of SMC and LMC stars. 
Further studies on a larger scale of the environment, and IR properties of the stars are needed. 
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1. Introduction 

We present the discovery and the study of two Small 
Magellanic Cloud blue irregular variable stars. We noticed 
that during EROS-2 SMC data base mining for microlens- 
ing by Palanque Delabrouille et al. (1998) 10 irregularly 
blue variable objects were detected (see their Sect. 3 to 
describe the method). Among them, two have the same 
variability characteristics as the Galactic Herbig Ae/Be 
stars (Herbig 1960, Bibo & The 1991, Waters & Waelkens 
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1998 and references therein), and were chosen for deeper 
analysis, similar to our previous studies of HAeBe candi- 
dates in the LMC (Beaulieu et al. 1996, and Lamers et 
al. 1999, de Wit et al. 2001a, BL, LBD, DBL hereafter). 
We discuss their photometric and spectroscopic proper- 
ties, and we present several arguments which might lead 
to the interpretation of these stars as pre-main sequence 
(PMS) HAeBe stars. If confirmed, these stars would be the 
first SMC pre-main sequence stars. Following the nomen- 
clature introduced in BL, we will call the stars EROS 
SMC HAeBe Candidates or ESHC stars. 

A large scale systematic search over the full EROS 2 
SMC database as done in the LMC by de Wit et al. 
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(2001a) will be presented in the near future (de Wit et 
al. 2001b). 

2. The observations 

2.1. Searching for blue irregular objects in the SMC 

The field has been observed as part of the ongoing EROS- 
2 microlensing survey. The set-up consists of a Im F/5 
Ritchey-Chretien telescope with two 2x4 CCD detectors 
in different focal planes covering a field of 0.7 x 1.4 de- 
gree. The light is split by a dichroic allowing simultaneous 
observing in two non standard bandpasses [Ve = 420-720 
nm and Re =620-920 nm). 

During the search for microlensing in the ER0S2 SMC 
data base (5.3 million stars coming from 10 square degrees 
of the densest part of the SMC) of one year observations 
(Palanque-Delabrouille et al., 1998), 10 blue irregular vari- 
able objects were discovered which are located in the HRD 
to the right of the main sequence. Two of these displayed 
irregular photometric variability with different time scales, 
superimposed on long term variability. The light curves 
resemble those of the HAeBe stars. 

The stars, called ESHCl and ESHC2, have the follow- 
ing coordinates : 

ESHCl: RA= 00:53:03.0 and DEC== -73:17:59, J2000 
ESHC2: RA= 00:52:32.6 and DEC= -73:17:08, J2000. 
The findings charts in V are given in Fig 1. 

The full Ve and Ve — Re light curves of theses two stars 
are shown in the two upper panels of Fig 2. The two lower 
panels show the low amplitude very irregular photometric 
behaviour in Ve and Re for both stars. The uncertainty 
in the photometry is typically cVe ~ cRe ~ 0.02 mag. 

The light curves of both stars show a long time scale 
behaviour (several hundred days AVe — 0.2 — 0.4 mag) 
with superimposed short time scale irregular variability 
(several days, AVe = 0.1 mag). They both have a bluer 
when fainter behaviour for the long time scale variation 
and the short time scale variation. The very irregular pho- 
tometric behaviour on short time scale is characteristic of 
Herbig AeBe stars, but may also be compatible with clas- 
sical Be stars. We will come back to this in Sect. 5. 

2.2. The photometry 

We noted that ESHCl is catalogued by Meysonnier & 
Azzopardi (1993) as an emission line star in their objec- 
tive prism survey, and identified as Lin 232 in the Simbad 
database. Despite the fact that ESHC2 was in the field 
observed by Meysonnier & Azzopardi, it was not detected 
as an Ha emitting object. 

On January 7, 1998 (JD 2450821.39) we have obtained 
BYRlHaRaoff photometry at ESQ La Silla with DFOSC 
on the Danish 1.5m for the SMC field, with exposure times 
of 60 sec, 30 sec, 30 sec, 30 sec, 600 sec and 600 sec respec- 
tively. The filters are the respectively numbered ESO 450, 
451, 452, 425, 693, 697. The CCD is a LORAL 2048x2048 
CCD with a pixel scale 0.39 arcsec. The data have been 



reduced using QUYLLURWASI, the pipeline developed 
around DoPhot for the PLANET collaboration. The ab- 
solute calibration has been done using Landolt standards 
taken during the observations, leading to absolute calibra- 
tion to 2 % for stars brighter than V=17 mag, and of the 
order of 5 % for stars of V=18 mag. 

Finally, using WCSTOOLS and the MACS catalogue 
(Tucholke et al., 1996), we performed astrometry of the 
reference frames and the subtracted frames. Comparison 
with other catalogues indicates that we have an accuracy 
better than 0.5 arcsec. 

We also searched for the program stars in the pub- 
licly available OGLE BVI catalogues for the SMC of 
2.2 million stars. These magnitudes are mean mag- 
nitudes over the observing period covered by OGLE 
(http://www.astrouw.edu.pl/ ^ ftp/ogle/). The photom- 



etry is summarized in Table 1. 

Moreover, narrow band imaging in Ha, Haoff was 
obtained with EF0SC2 mounted on the ESO 3.6m at 
La Silla on August 18, 1999 (Juhan day 2451408.79) for 
the field containing the targets with exposure times of 
600 sec. Flatfielding, bias subtraction, cosmic subtraction 
were performed with MIDAS 98NOV in a standard way. 
Although the seeing was poor (1.3 arcsec), these images 
are useful to search for Ha emitters in the field. We built 
subtracted images (Ha — Haoff) using ISIS, the new im- 
age subtraction package with non constant kernel (Alard 
2000). The procedure is the following : we choose the Ha 
image as being the template. We re-map the Haoff image 
by bicubic spline interpolation from the Ha image. Then 
we adjust the parameters of ISIS to our images, and do 
several tests with different degrees to fit the variation of 
the convolution kernel over the image. We also performed 
several tests to fit the background variations. The best re- 
sults are obtained with a degree 1 variation in both the 
kernel and the background. 

We present images of the (Ha — Haoff) /V Ha flux in 
the lower left panel and the right panel of Fig. 1. It gives an 
idea of the noise statistics compared to the photon noise 
: the gray scale has been adjusted so that black means 
4 sigma above photon noise level. It is clear that both 
ESHC objects are strong Ha emitters. They lie on a faint 
HII region illuminated by nearby hot stars. Notice that 
at 5 arcsec (1.5 pc on the sky) to the North West from 
ESHC2, there is another Ha emitter. Moreover, there are 
several Ha emitters clustered in the upper left corner of 
the field. This indicates that we are most likely in an area 
where star formation was active recently. 

2.3. The spectroscopic observations 

Spectroscopic observations were obtained with EF0SC2 
mounted on the ESO 3.6m at La Silla on August 18, 1999. 
The gratings nr. 7 and nr. 9 were used to cover the full 
spectral range from 3200 to 7000 A. We used a 1.2 arcsec 
slit, and we binned the data to a resolution of 6 A. Multiple 
exposures of the two stars were made with the grating nr. 
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Fig. 1. The upper left panel is the finding chart (3.25 arcmin x 1.65 arcmin) in the Haoff band of the two targets 
ESHCl and ESHC2, labeled by the left and the right cross respectively. North is up, East is left. The two other are 
{Ha — Haoff)/ Vila images. These figures show only the Ha point sources and the Ha emitting regions (stripes are 
related to data reduction problems like bad columns or saturated stars). The lower left panel is the area shown on the 
finding chart, the right panel gives the full field observed. One can notice the presence of an HII region, and clusters 
of Ha emitters. 



Name 


Vd 


(B-V)d 


(V-R)d 


(V-I)d 


Ha — Haoff 


Vog 


(B-V)og 


(V-I)og 


ESHCl 


15.09 ± 0.02 


0.03 ± 0.02 


0.06 ± 0.02 


0.11 ±0.02 


-0.42 ± 0.03 


15.07 ±0.01 


-0.00 ± 0.03 


0.13 ±0.04 


ESHC2 


17.03 ± 0.02 


-0.03 ± 0.02 


0.05 ± 0.02 


0.10 ±0.02 


-0.58 ± 0.03 


17.01 ±0.02 


-0.05 ± 0.03 


0.07 ± 0.04 



Table 1. Photometry obtained for ESHCl and ESHC2 with DFOSC on the Danish telescope at ESO La Silla on 
January 7, 1998 (indicated by a "d") and from the published catalogue of average SMC BVI photometry by OGLE 
(indicated by a "og"). 



7 (1200 sec exposure), and grating nr. 9 (900 sec exposure) 
for each star. The spectra were reduced in a standard way 
using MIDAS 98NOV release and the LONG context. A 
signal to noise ratio of 200 was obtained for ESHCl and 50 
for ESHC2 in the blue part of the spectrum. A somewhat 
lower signal to noise ratio was obtained in the red part of 
the spectrum: 100 for ESHCl and 30 for ESHC2. 

3. Discussion of the spectra 

The full spectra are presented in the Fig. 3 with ESHCl 
in the upper panel and ESHC2 in the lower panel. The 
region between 3800 < A < 4500 A is enlarged in Fig. 

4. The H/9 and Ha profiles are shown in Figs. 5 and 6. 
We immediately notice that these two stars have stellar 
spectra with small Balmer jumps, clear Balmcr lines in 
absorption and strong Ha emission. This suggests spectral 
type Be. 

In order to improve the estimated spectral type, we 
decided to adopt the same approach as in LED. We first 



looked for the presence or absence of classification lines, 
Balmer lines, Hel, Hell, CHI, OH, Mgll, Sill, SilV. When 
it was possible, we measured the position and equivalent 
width of the lines. It is clear that the metal lines are 
weak, as expected for metal poor stars in the SMC, and 
marginally detected if detected at all. Therefore, in the 
following we will concentrate on the Balmer lines and the 
Helium lines. The results are reported in Table 2. In the 
next two subsections we will discuss the spectral type of 
the two program stars. 

3.1. The spectrum of ESHCl 

ESHCl shows Ha in emission, H/3 is partly filled with 
emission, all other Balmer lines from H7 to H12 and 
Helium I (3819, 3888, 4026, 4387, 4471, 4713, 4921, 5015, 
6678 A) are present, and clearly seen in absorption (Fig. 
4). This star has a clear Balmer jump. The Hell line at 
4686 A is absent. There may be a weak contribution of 
the CHI/OH blend at 4650, but if present, its equivalent 
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Fig. 2. Light curves of the two targets. The two upper panels shows the full Ve and Ve — Re light curves for ESHCl 
and ESHC2. The two lower panels show zooms on a part of the light curves. Notice the low amplitude irregular 
photometric variability. The time is given as JD-2447892.5. A typical error bar is given in the lower right of the two 
upper panels. 



width would be less than 0.1 A. There may be a trace of 
on at 4662 A. There may also be a trace of SilV at 4116, 
but the SilV counterpart at 4089 A is in the wing of He. 
Therefore the presence of these lines is inconclusive. The 
Mgll doublet at 4481 A is not detected. 

Wc compare the data with the tables from Didelon 
(1982) of the equivalent width as a function of spectral 
type for several lines. Most of the metal lines do not give 
strong constraints in this particular case. However, the 
presence and strength of the H and Hel lines suggest that 
this star is a Be star with a luminosity class III, IV or V. 
The spectral type and luminosity class is best constrained 
by the Hel 4387 line, which indicates a luminosity class III, 
and a spectral type B4-6. The Hcl 4471 line also implies 
that the star has a type between B4 and B6, but Hcl 4026 
line implies a slightly earlier spectral range of Bl to B4. 
Based on these lines, we classify this star as a B4e HI. 

3.2. The spectrum of ESHC2 

ESHC2 shows Ha and H/3 in emission. All other Balmer 

lines, from H7 to H13, and Hel (3888, 4026, 4143, 4392, 
4471 A) are present and clearly in absorption. This star 



has a weaker Balmer Jump than ESHCl. The Hell line at 
4686 A and the Mgll doublet at 4481 A are absent. No 
features are seen at 4650 A from the CIII/OII blend. Just 
Uke for ESHCl, we used the tables of Didelon (1982) on 
the Hel lines at 4471, 4026, 4387 and 4120 A to determine 
the spectral type and luminosity. Based on the compari- 
son of the observed line strength with Didelon's tables we 
classify this star as a B2e IV-V. 



3.3. Spectral type from the Q method 

We can also estimate the spectral type from the pho- 
tometry, using the relation between the reddening-free 
Q-parameter and the spectral type. For this purpose we 
adopt the transmission curves of the U, B and V filters and 
convolve these with the observed spectra to derive the U , 
B and V magnitude in the Johnson system. The resulting 
colours are for ESHCl [U - B = -0.45, B -V = 0.04) 
and for ESHC2 [U - B = -0.62 and B - F = 0.20). 
We apply the seminal Q method (Johnson & Morgan 
1953) to determine reddening-free colour Q = {U — B) — 
0.72(B - V). Wc derived Q = -0.47 for ESHCl and 
Q = —0.76 for ESHC2. Using the relation between Q and 



Beaulicu ct al.: Small Magellanic Cloud blue irregular variables 



5 



Star ESHCl 


Star ESHC2 




position 


FWHM 


EW 


Flux 


position 


FWHM 


EW 


Flux 


Identification 


(A) 


(A) 


(A) 




(A) 


(A) 


(A) 






3822.8 


6.0 


-0.54 


-17 










Hel 3819 


3892.0 


10.4 


-3.2 


-106 


3889.8 


18.0 


-2.7 


-16 


Hel 3888 


4012.2 


6. 


-0.4 


-12 










Hel 4009 


4029.4 


7. 


-0.8 


-23 


4028.4 


11 


-0.8 


-4.5 


Hel 4026 


4146.5 


7.6 


-0.6 


-17.5 


4147.7 


7 


-0.7 


-4. 


Hel 4146 


4391.5 


9.3 


-0.74 


-20.2 


4392.0 


7.6 


-0.8-1.0 




Hel 4387 


4475.1 


8.8 


-0.83 


-22 


4473.8 


6 


-0.6 


-3 


Hel 4471 


4717.3 


7.5 


-0.28 


-6.8 










Hel 4713 


4925.7 


6.8 


-0.62 


-13 


4925.0 


7.4 


-0.7 


-3.4 


Hel 4921 


5019.6 


6.2 


-0.23 


-4.65 










Hel 5015 


6684.3 


10.9 


-0.6 


-5.6 










Hel 6678 


5894.5 


12.6 


-1.0 


-13 


5890.6 


14 


-1.0 


-3.2 


Na 


3737.3 


6.0 


-1.1 


-33 










Hw 3734 


3753.0 


7.5 


-1.7 


-51 










Hu 3750 


3773.4 


7.9 


-2.1 


-65 


3773.1 


13 


-2 


-11 


Hii 3770 


3800.9 


9.12 


-2.58 


-81 


3800 


14 


-2.3 


-13 


J^io 3798 


3838.4 


8.8 


-2.5 


-81 










Hel 3838 + Hg 3835 


3936.8 


asym 






3934.2 


9.7 


-0.9 


-5.4 


Ca H 


3972.7 


10.4 


-3.4 


-110 


3972.1 


13.9 


-4.0 


-24 


He / He K 


4105.1 


9.7 


-2.33 


-70 


4105.4 


14.1 


-2.7 


-16 


m 


4344.2 


12.7 


-2.54 


-69 




filled 








4865.5 


partly filled 


4864.5 


9 


2 


9 


H/3 


6567.8 


7.5 


6.2 


58 


6567.4 


10.1 


18.5 


40.4 


Ha 



Table 2. Lines identified in ESHCl and ESHC2. We give the measured wavelength of the lines (columns 1,5), the 
Full Width at Half Maximum (given in columns 2, 6), the equivalent width (column 3, 7) and the measured Flux (in 
columns 4, 8). 



spectral type for Be stars by Halbedel (1993) we derive 
a spectral type B5 and B2 respectively. This agrees very 
well with the spectral types of B4eni and B2eIV-V de- 
rived from the H and Hel lines above. 

A word of caution may be necessary here. The Q ver- 
sus spectral type relation that we adopted is derived for 
Galactic stars. The metal poor SMC stars, which have 
less line blanketing and line blocking may have different 
colours than Galactic stars of the same T^g, and therefore 
different temperatures than the Galactic stars of the same 
Q-value. However, comparing the colours and magnitudes 
of the Kurucz model atmospheres (1979) for Z = 0.02 
(Galactic) and Z = 0.002 (lower than SMC), we find that 
the difference is very small. Given the uncertainty in the 
method and the fact that the H and Hel lines indicate the 
same type as the Q-method, we will adopt these spectral 
types and their corresponding values of Teg from Galactic 
stars. 

4. Determination of the stellar parameters T^g, L, 
M, R 

4.1. ESHCl 

We adopt a spectral type B4 HI for ESHCl as the most 
reasonable estimate. From the OGLE photometry the 
mean magnitude and colours are V — 15.08 ± 0.01 and 
{B -V) = 0.00 ± 0.03. The intrinsic colour of a B4III 



star is {B - V)o = -0.19 (Lang 1991), which results 
in E{B -V) = 0.19 ± 0.03. We make the assumption 
that the colour difference is due to obscuration by nor- 
mal dust with Rv = 3.3 (Cardelli et al., 1989). This 
results in an intrinsic magnitude of Vq = 14.45 ± 0.11. 
Adopting a distance modulus of 18.94 ± 0.04 (Laney & 
Stobie 1994) to the SMC, we find My = -4.49 ± 0.15. 
With a bolometric correction of -1.45 for spectral type 
B4in (Lang 1991) this yields = -5.94 ± 0.15 or 

log(_L/LQ) = 4.23 ± 0.06. The effective temperature of 
a B4III star is log T^s = 4.20 ± 0.02 (Lang 1991) and 
so log(i?/i?Q) = 1.24 ± 0.03 and R = 17 A ± 1.2 i?©. 
Following Lang, we use an empirical Mass Luminosity re- 
lation log(M/MQ) = 0.48 - 0.105Mboi and derive a mass 
of M pa 13M0. Similar results are obtained when we de- 
rive the parameters with the photometry we obtained from 
ESQ telescope. 

4.2. ESHC2 

We adopt a spectral type B2 IV-V for ESHC2 as the 
most reasonable estimate. From the OGLE photometry 
the mean magnitude and colours are V — 17.01 ± 0.02 
and {B — V) = —0.06 ± 0.03. Adopting the intrinsic 
colour {B — V)o = —0.24, an effective temperature of 
log Tcff = 4.34 ± 0.02 and a bolometric correction of 
BC = —2.35 from Lang (1991) we find the following stel- 
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spectrum of the two targets taken with Fig. 5. The Hl3 line and Hel lines at 4713 A and 4925 A 
the ESO 3.6m. The spectra have a resolu- of ESHCl and ESHC2. 



lar parameters: Vq = 16.42 ± 0.12, My = -2.52 ± 0.16, 

Mboi = -4.87 ± 0.16, log(i/L0) = 3.84 ± 0.06 and 
i? = 6.1 ± IRq. We get a mass of M lOM©. Again, 
similar results are obtained when we derive the parameters 




He, 



Hd 



H7 



3800 



4000 4200 
wavelength (A) 



4400 



P 




6500 



6550 6600 6650 

WAVELENGTH (A) 



6700 



Fig. 6. The Ha line of ESHCl and ESHC2. 

with the photometry we obtained from the ESO telescope. 
The results are summarized in Table 3. 

5. Discussion 



Fig. 4. The blue part of the spectrum of ESHCl and We have presented two stars with the following character- 

ESHC2. Balmer lines He, RS, H7 are labeled, and the istics: 

solid lines vertical indicate Hel lines at 3822, 3892, 4012, (1) early B spectral type with strong Ha emission, 

4026, 4146, 4387 and 4471 A. ^2) irregular photometric behaviour with time scales of 
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Name 


Type 


logTeff 


log(L/L0) 


{R/Rq) 


ESHCl 


B4 III 


4.20 ± 0.02 


4.23 ± 0.06 


17.4 ±1.2 


ESHC2 


B2 IV-V 


4.34 ± 0.06 


3.84 ± 0.06 


6.1 ±1 



Table 3. Summary of stellar parameters for ESHCl and 
ESHC2. 



20-200 days, 

(3) a bluer photometric colour, during brightness mini- 
mum, 

(4) located within a projected distance of 20 pc of each 
other in a weak HII region, 

(5) positioned in the HR-diagram to the right of the Main 
Sequence (see below). 

The five points suggest that the stars might be Herbig 
AeBe or classical Be stars. If they are PMS stars, then 
they are the first discovered Pre-Main Sequence stars in 
the SMC (see also Hutchings & Thompson 1988). 

We rule out that the objects arc B[e] supergiants, as 
no forbidden line emission was detected. Moreover the 
ESHCs have a lower luminosity class, at variance with the 
supergiant nature of B[e] stars in the Magellanic Clouds 
(Lamers et al. 1998). 

5.1. Comparing the properties of ESHCl and ESHC2 
with classical Be and HAeBe stars 

To make a distinction based on visual properties between 
the two possible classification is precarious as many prop- 
erties are common to both classical Be stars and Herbig 
Be stars. However we might find some indications, by hav- 
ing a closer look at the characteristics mentioned above. 
The known physical processes (ie. bound- free and free- free 
continuum emission, noted bf and ff respectively hereafter) 
which occur in the CS environment of classical Be stars 
allow predictions with respect to the Ha emission and 
colour excess E{B — V). In contrast, the activity exhibited 
by Galactic HAeBe stars remains fuelling active debate, 
and therefore no quantitative predictions can be made. 
- The Ha emission of Classical Be stars is correlated 
with the observed colour excess (Dachs 1988). It is ex- 
pected that a larger contribution of the bf-ff excess emis- 
sion will generate larger Ha emission accompanied by 
a larger colour excess. Early type Galactic classical Be 
stars have a maximal colour excesses of E{B — F) ~ 0.10 
(Waters et al. 1987). The derived E{B - V) for ESHCl 
and ESHC2 are 0.19 and 0.18. If we correct for the average 
extinction towards the SMC {E{B-V) = 0.09, Massey et 
al. 1995), we find that the colour excesses of ESHCl and 
ESHC2 are similar to the maximum observed for the early 
type Galactic Be stars (Waters et al., 1987). However, 
comparing the colour excess with the Ha emission for clas- 
sical Be stars done by Dachs (1988) we noticed that the 
colour excess is much larger than expected on the basis of 
the Ha emission. This might be attributed to additional 
circumstellar extinction. In general Galactic HAeBe have 



colour excess which exceeds the level observed for classical 
Be stars. 

- The irregular brightness variability of Classical 
Be stars has been shown to be largest among spectral 
types B1-B2, on the basis of the four years (similar to 

the ER0S2 data) of observations of Hipparcos (Hubert & 
Floquet 1998). Amplitudes of variability in the Hipparcos 
passband were measured to have maximum values of 
AHp = 0.3, associated with outbursts on time spans 
of hundreds of days. This variability is accompanied by 
small amplitude (< 0.1"*) variation on time scales of 
days, often periodic (e.g. Pavlovski et al. 1997, Hubert 
& Floquet 1998). The amplitude of the variabihty of the 
ESHCs is mildly larger the classical Be stars, however 
the Hipparcos filter is wider and extends further to the 
blue part of the spectrum than the ER0S2 Ve passband. 
Bf-ff emission has its largest contribution in the red 
part. Therefore the amplitude in the ER0S2 passband 
is expected to be mildly larger than the Hipparcos 
passband. The superposed short time scale variability 
of the ESHCs (^ 20 days), is however not reported for 
Be stars by Hubert & Floquet (1998). This time-scale of 
variability is however observed among the HAeBe stars 
(Herbst & Shevchenko 1999). 

- The observed colour variability of the ESHCs 
can be interpreted in terms of physical processes, which 
are thought to occur in the CS environment of classical 
Be stars (Dachs 1982, Waters et al. 1987). However 
this explanation is not unique. Just like the two SMC 
HAeBe candidates, most of the LMC Pre-Main Sequence 
candidates also show bluer-when-fainter behaviour. In 
LBD we discuss the possibility of explaining the observed 
photometric and colour variability in terms of dust 
obscuration and the existence of a reflection nebula. The 
existence of a reflection nebula is one of the original 
selection criteria of Herbig (1960). The total flux of 
HAcBc stars at large distances (like in the LMC and 
in the SMC), will be contaminated by the unresolved, 
blue reflection nebula (typical sizes of Galactic reflection 
nebulae ~ 0.25 pc will be less than 1 arcsecond in the 
Magellanic Clouds). The nebula will also cause the colour 
excess to be small. In LBD it was concluded that one 
would preferentially observe Magellanic Cloud HAeBes 
having a bluer-when-fainter colour behaviour due to dust 
obscuration. 



- The close proximity of the two ESHCs to each other 

would mildly support the HAcBc classification. However, 
we cannot exclude the chance nearby projection of two 
irregularly variable classical Be stars. The H II region does 
indicate the presence of hot stars and thus a region which 
might have recently witnessed star formation. 
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5.2. Position in tlie HR-diagram 

We have seen in the previous subsection that the (visual) 
evidence does not point towards one particular classifica- 
tion (classical Be or Herbig Be) of the ESHCs. In the two 
panels of Fig. 7 we present an HR-diagram. The purpose 
is to compare the derived stellar parameters of the ESHCs 
with the Classical Be stars and the Herbig Ae/Be stars. 
The ESHCs are presented as filled symbols. The Galactic 
HAeBes are plotted as crosses, and the open symbols are 
their counterpart candidates in the LMC, the ELHCs. The 
full line is the Main-Sequence. 

The dotted lines in the upper panel of Fig. 7 repre- 
sent the "birthline" for pre-Main Sequence stars for two 
different proto-stcUar mass accretion rates. The lower dot- 
ted line is the canonical Galactic birthline for a accretion 
rate of lO~^M0yr~^, and higher dotted line is the birth- 
line for an accretion rate 10~'^Mq yT~^ (Palla & Stabler 
1993). The crosses are Galactic HAeBe stars adopted from 
Berrilli et al. (1992). The observed upper limit for Galactic 
HAeBes agrees approximately with an accretion rate of 
10~^Mq yr~^. We see, as noticed in BL and LBD, that 
some of the LMC HAeBe candidates (open symbols) are 
more luminous by a factor 3 to 10 than the Galactic birth- 
line. In comparison the luminosity of SMC star ESHCl 
complies with the high luminosity found for the LMC 
HAeBe candidates. 

The luminosity of the birthline depends on the accre- 
tion rate. We found that some of the ELHCs and ESHCl 
agree better with the birthline corresponding to a fast 
accretion rate of 1O~^M0 yr~^. Just based on these two 
stars, we cannot compare birthlines in LMC and SMC yet, 
but it is clear that these two pre- main sequence candidates 
in the SMC are above the Galactic birthline, and are very 
similar to the LMC pre-main sequence candidates. If in- 
deed they are pre-main sequence stars, it would support 
our previous conclusion that, at low metallicity, we can 
observe more massive stars in pre-main sequence phases 
than in our Galaxy. 

Palla (2000, private communication) has argued that 
massive Galactic PMS may also be above the birthline for 
an accretion rate of lO~^M0yr~-'^. If that is the case, the 
high location of ELHCs and ESHCs in the HRD may not 
be necessarily related to the low metallicity. On the other 
hand the massive PMS stars noted by Palla are closer to 
the Main Sequence than the ELHCs and ESHCs. So, in 
the temperature range 10000 < Tg// < 15000, the LMC 
and SMC HAeBe candidates are definitively above the 
Galactic HAeBe stars. 

On the other hand, we consider the case that they can 
be classical Be stars. The dashed lines in the lower panel of 
Fig. 7 indicate the upper limit for main sequence and giant 
classical Be stars. These two upper limits have been de- 
rived from Zorec & Briot (1997). In this work the average 
V-Magnitude excess of classical Be stars has been derived 
as function of B spectral subtype. We added this excess to 
the the absolute visual magnitude (also given in Zorec & 
Briot 1997), in which we included the dispersion in 



order to obtain the upper limit. The effective temperature 
and bolometric correction were adopted from Schmidt- 
Kaler (1982). The figure shows that the stellar parame- 
ters of ESHC2 are compatible with the derived classical 
Be upper limits for main sequence stars, and ESHCl with 
the upper limits for giants. 

It is not possible to unambiguously distinguish be- 
tween HAeBe and classical Be interpretation for the two 
stars! 

6. Conclusions 

We have identified two blue irregular variables in the 
SMC with spectral type B4IIIe (ESHCl) and B2IV-Ve 
(ESHC2), with strong Ha emission. These stars could 
be either classical Be stars or pre-main Sequence HAeBe 
stars. We find evidence for both interpretations; 

— The Ha emission of the ESHCs stars is much larger 
than that of the classical Be stars with the same colour 
excess. 

— The short timescale variability of the ESHC stars of 
about 20 days is similar to that of HAeBe stars, but 
not found in classical Be stars. 

— The bluer- when-fainter colour variability could be ex- 
plained by variable free-free and bound-free emission 
(as in classical Be stars) or by variable extinction of a 
star with a blue scattering nebula that is not resolved 
(as expected for HAeBe stars at the SMC distance). 

— The presence of a nearby H II region and the proximity 
of clusters of iJa-emitters indicates recent star forma- 
tion, which would favour the HAeBe interpretation. 

So we find a mild preference for the interpretation of 
the ESHC stars as HAeBe stars. On the other hand, the 
possibility that the stars are classical Be stars with abnor- 
mal photometric properies cannot be excluded. Infrared 
observations and high resolution observations are needed 
to confirm or deny one of the two possibilities, because 
HAeBe stars are known to have an IR excess due to dust 
emission and HAeBe stars are often surrounden by a re- 
flection nebula. 

If the stars turn out to be HAeBe stars indeed, they are 
the first pre-main sequence stars found in the SMC. Their 
magnitude then indicates that they are located above the 
Galactic birthline in the HR-diagram, which suggests a 
faster accretion rate than their Galactic counterparts. 
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log(Teff) 

Fig. 7. Both panels are HR diagram for Galactic (crosses, 
Berrilli et aL 1992), LMC (open circles, Lamers et aL 1999) 
Herbig Ae Be candidates and the two targets of this study. 
The main sequence is indicated as a line. The upper panel 
shows also the two birthlines from Palla & Stabler (1993) 
for mass accretion rates of 10~^Mq yr~^ (long dashed 
lines) and 1O~'*M0 yr~^ (short dashed lines). The lower 
panel shows (long dashed lines) the upper limit for main 
sequence Be stars and giant Be stars, adopted from Zorec 
& Briot (1997). 



